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1 Introduction 



Understanding the nature of black hole thermodynamics may teach us lessons about 
quantum gravity. The underlying phase structure can be useful not only in this regard but 
for other purposes as well. For example, with the advent of AdS/CFT correspondence, the 
known phase structure of large black holes in asymptotically Anti-de Sitter (AdS) space 
can be used to understand various physical phenomena in other branches of physics. An 
illustration of this is that the Hawking-Page transition for AdS black hole 'evaporating' 
into regular "hot empty AdS space" at certain temperature [1] can enhance understanding 
of the confinement- deconfinement phase transition in large N gauge theory [2] . 

A large part of the phase structure of an AdS black hole [HI H] is actually not unique to 
the black hole in asymptotically AdS space but shared universally by suitably stabilized 
black holes/branes, say, in asymptotically flat space O El [TJ |8] , even in the presence of a 
charge q. For example, a chargeless (suitably stabilized) asymptotically fiat black p-brane 
can also undergo a Hawking-Page transition at certain temperature, now evaporating into 
a regular 'hot fiat space' instead. When q ^ 0, there exists also a critical charge qc and 
for q < qc, the phase diagram universally contains a van der Waals-Maxwell liquid-gas 
type phase structure along with a first-order phase transition line ending at a second- 
order critical point with a universal exponent for the specific heat as —2/3 when q = q^.. 
This universal phase structure may hint holography even in asymptotically fiat space, as 
pointed out in [5]. 

However, unlike an AdS black hole, an isolated asymptotically flat black hole/brane 
is unstable due to its Hawking radiation and needs to be stabilized first before one can 
discuss the equilibrium thermodynamics. To establish its stability, the standard practice 
is to place such a system inside a finite spherical cavity [9] with its surface temperature 
fixed. In other words, a thermodynamical ensemble is considered which can be either 
canonical or grand canonical, depending on whether the charge inside the cavity or the 
potential at the surface of the cavity is fixed [10] . In this paper, our focus is the canonical 
ensemble, i.e., the charge inside the cavity is fixed, and in particular our main interest is 
to study the phase structure when the fiux/charge inside the cavity is fixed but non-zero. 

When the ensemble temperature drops below a certain minimum value, there appears 
a puzzle of missing phases, in certain CcLSGS, clS it is not clear where the system would be 
in the absence of globally stable phases (No such issues arise in grand canonical ensemble, 
however, [3 [8].). Furthermore, it is not known whether there exists a new globally stable 
phase in the present setup other than what have been discussed in the literature so far 
[3 ini [7] . This was noticed in [5l [6] for the charged black hole and in [^ for the charged 
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black p-brane when d <2. 

In this paper, we will resolve the above puzzle and give the new phase(s) for the D- 
dimensional asymptotically fiat stabilized compact black p-brane (for p > 0) □ by finding 
two missing phases, namely, the regular 'hot bubble', due to the existence of "bubble of 
nothing" [11], and the extremal brane, each carrying the same fiux/charge as the black p- 
brane. The bubble or the extremal brane [121 ESI IHl ES] each can have an arbitrary period 
/3 in Euclidean time, analogous to the 'hot fiat space' in the chargeless case. As such the 
black p-brane can make a transition to this bubble or the extremal brane, depending on 
whose free energy is smaller, giving an analog of Hawking-Page transition. As a result, the 
underlying phase structure is greatly enriched and many new phase transitions between 
black branes, bubbles and extremal brane are revealed, giving a rather complete phase 
structure in this setup. 

This paper is organized as follows. In section 2, we will present the basic setup for 
phase structure of black branes which will be discussed in the following sections. In 
section 3, we discuss the phase structure for the special zero fiux/charge case as a warm 
up exercise. Here we see how the inclusion of bubble phase will enrich the previously 
known phase structure. Section 4 is the main focus of this paper and we will resolve the 
aforementioned puzzle for d < 2 and find a new global stable (bubble) phase for d > 2, 
therefore giving a necessary completion of underlying phase structure and various new 
phase transitions including the analog of Hawking-Page transition among other things. 
We discuss the results obtained in this paper and conclude in section 5. 

2 The basic setup 

For the purpose of this paper, let us consider the D-dimensional black p-brane metric in 
Euclidean signature as [161 E] , 

dsl, = A+AZ^df + Af^ {dx^f + A?^ J^'^dx'y + A;i A^^'dp' + p'A^dfiJ^^ (1) 

1=2 

where A± = 1 — {r±/pY, with r±. (r^. > r„) related to the mass and the charge of the 
black p-brane. The horizon occurs at p = r+ while the curvature singularity at p = r_. 

■^For p = 0, there is no bubble phase. The only phases are non-extremal brane and extremal one, and 
the phase relation between the two for p > discussed in the text will hold true also for this case. For 
simplicity, we focus from now on only for p > case. 
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Here d=l+p, d = D — d — 2 and 'a' is the dilaton coupling defined by 

for supergravity with maximal supersymmetry. To have a large but finite Euclidean action 
for the black brane, the brane directions x*, with i = 1, . . . ,p should be compact (In the 
metric, coordinate is explicitly isolated for the purpose of constructing bubble solution 
later). For the metric ([1]) to be free from conical singularity at p = r+, 't' coordinate must 
be periodic with periodicity 



47rr_i. / r 



d 



the inverse of the black brane temperature at p = 00. With this, the inverse of the local 
temperature at p is 



. (4) 



However, for the corresponding extremal brane, 't' coordinate can have an arbitrary period 
[T2| [T3| [T5] and this is crucial for the phase transition we will discuss later. Here the 
physical radius p = A" ^^'^p can be read from the metric ([1]) and so the physical parameters 
f± = A"L^^'^r±. On the other hand, the coordinate x^, like other compact coordinates, 
has arbitrary local periodicity. By renaming the coordinates t and t — )■ — in 

the black p-brane configuration given above, we can obtain the bubble carrying the same 
flux/charge in Euclidean signature with its metric given a^ 

dsl^ = A^df + A+ Al^ {dxy + A?^ Y.'^dxy + A;^A^"'rfp2 + p^Apdil}^^, (5) 



^Though this provides a simple means to obtain a bubble solution, it can also be obtained by other 
means, for example, by solving the underlying equations of motion. In other words, once the solution 
is obtained, we can forget about its connection with the original black brane. Even in the present 
context of considering possible allowed phases with the same boundary data, the parameters r_)_ and r_ , 
characterizing the solution, can be different in both the bubble solution jSj and the black brane ([T]). As 
such, the period of (or t) in the bubble case is not necessarily the same as that of t (or x^) in the black 
brane case. Only when the r+ and r_ are set to be the same in both cases, the two periods will be the 
same which is just a special case. This recognition, as discussed later in the text, is crucial for the role 
of the bubble phase in the phase diagram uncovered in this paper. 
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which is defined only for p > r+ and is regular if has a period of 

at p = r4.. However, in this case the periodicity of the coordinate 't' as well as the other 
spatial compact coordinates remain arbitrary and so it can be in thermal equilibrium with 
the cavity in any temperature. 

Now in order to study equilibrium thermodynamics [18] in canonical ensemble, the 
allowed configuration (black brane or extremal brane or bubble or coexistence of two 
or more of them) should be placed in a cavity [9] with fixed physical radius pB- The 
other fixed quantities are the cavity temperature the physical periodicity of x^, 

i.e., R (also the physical sizes of other compact directions), the dilaton value and the 
charge/fiux enclosed in the cavity Qd- In equilibrium, these fixed values are set equal to 
the corresponding ones of the allowed configuration enclosed in the cavity. For example, 
we set the charge 

where * denotes the Hodge duality and -F[p+2] = id{{f^f^Y/p'^~^^)dpAdtAdx^ A ■ ■ ■ Adx^, 
the field strength for the configuration considered. In the charge expression Qn denotes 
the volume of a unit n-sphere, k is a constant with 1/(2k^) appearing in front of the 
Hilbert-Einstein action in canonical frame but containing no asymptotic string coupling 
Qs- We have also set = e''^^^^^ = QgA'^'^. In canonical ensemble, it is the Helmholtz 
free energy which determines the stability of the equilibrium states and is related to 
the Euclidean action by F = J^;//? in the leading order approximation. So, in order to 
understand the phase structure we will evaluate the action for the black (non-extremal) 
brane, the extremal brane and the bubble, with the above mentioned boundary data. 
Note that specifying the boundary data does not necessarily mean that they have the 
same f+. For example, in the extremal case, f+ = f_ are completely fixed by the given 
charge. 

The Euclidean action for the black p-brane has been evaluated in [7] by using standard 
technique and is given as, 

'A 
A 



.bi _ (^RVp-i^d+i -d 

— o, 2 



1/2 

(2 + (i)(=i) + J(A_A+)i/2_2(J. 



rf+l -d+l A" 
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with A± taking their respective value at p = pB- Since the Helmholtz free energy is given 
as Fbi = -Ebi — ^'S'bi, so we have J^' = /3i?bi — •S'bi, where i?bi is the internal energy and is 
the entropy of the black p-brane. Thus we identify the internal energy of the black brane 
on dividing the first term in (JHl) by /3 and can be checked to match the ADM mass per 
unit volume of the black brane [19] as — ^ oo. The second term in ([8]) is the entropy of 
the black brane. Note that we have written the usual compact brane volume Vp = RVp^i. 
The action for the bubble can be obtained simply from ([8]) by making the change 
13 ^ R and is given as 



rbb 




a/2 



(9) 



As for the bubble since there is no entropy, we have = /5£'bb and so we can identify the 
internal energy of the bubble on dividing its action by /3 and again can be checked to have 
the correct ADM mass per unit volume of the bubble as pB oo. However, for the black 
brane in canonical ensemble the local temperature 1/P{Pb) is determined by r^,Q(i,pB 
as given before and the periodicity of is arbitrary whereas for the bubble the local 
periodicity of x^, i.e. R{Pb), is similarly determined by the corresponding quantities and 
the temperature is now arbitrary. On-shell, they are all set equal to the corresponding 
fixed boundary values. 

Now for convenience, we will work instead with the reduced action defined as i^'^^{z) = 



Gq{z) with Gq{z) defined as. 



-bR 



{d + 2) 




z) 



1/2 



1/2 



(10) 



where we have defined z = [t^/PbY < 1 £^nd z = x for the black brane and z = y for the 
bubble since the two need not have the same f_|_. Further we have defined b = ^/{AttPb), 
R = R/ (AtiPb) and T = _R for the black brane and T = 6 for the bubble. Also, q = {Q*^/ psY, 
where = [(A/2K(5d)/(f^j+irf)]^^'^. In terms of these new fixed parameters we have. 



A, 



A_ = 1 - q^z, 1-^:^ = 1- q^/z^ 



fill 
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and they have been used in (fTOl) . In writing A_ we have used the fact that since Qd is 
fixed, f_ is not an independent parameter but can be expressed in terms of f^_. Given 
{Q*df'/^'+ = f_/f+ < 1, SO z > g. For extremal branes, f_ = f_|_ and so z = q. The 
reduced action for extremal branes is now 

^extremal _ ^^J^^ ^^2) 

determined completely by the boundary data. However, for non-extremal brane or bubble, 
we have a variable z lying between q < z <1. Now dGq{z) / dz\^^- = gives the equation 
of state m = 'mq{z), where, 

mJz) = -, — (13) 

If 

(fGg{z)/dz'^\^^. oc - dmg{z)/dz\^^- > 0, (14) 

ai z = z, where z is determined from m = mq{z), we get a local minimum of free energy. 
So the negative slope of mq{z) determines the local stability of the underlying system. In 
the above rh = b, mq{z = x) = bq{x) for the black brane and rh = R, mq{z = y) = Rq{y) 
for the bubble (Note that we always have fhl = bR.). 

The phase structure of Gq{z) corresponding to the black p-brane has been analyzed in 
[7j. The bubble also has exactly similar phase structure, however, the relevant quantities 
here are R and Rq{y) instead. For now, we just need to compare the free energies among 
the black brane, the bubble, at their respective global minimum, and the extremal brane, 
all with the same boundary data. The analog of Hawking-Page transition in either case 
and the final stable state will be determined by the smallest free energy of these phases. 
For this purpose we need the respective on-shell reduced free energy explicitly. It is 

f bl, bb / -\ 

F-^^'^^^) = -^-r^ = -RFqizl (15) 



with 



1/2 „2\l/2 



Fqiz) = 2\^\ +d\Y^j +dil-zf'(^l-^-j -2(rf+l). (16) 

In the above, we have used the on-shell condition rh = mq{z) with z lying between 
q < z < 1. Given the boundary data b and R, the on-shell free energy of non-extremal 
brane is F"(x) = —RFq{x), with X determined hj b = bq{x), while the on-shell free energy 



2\ 1/2 
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of bubble is F^^{y) = —RFq{y), but now with y determined by _R = Rq{y). So F^^{x) 
and F^^{y) actually have different dependence on their respective on-shell variable even 
though both have the same functional form —RFq{z) in appearance. So their profiles 
are in general different. For example, F^^{x = q) = Rdq (giving the extremal brane free 
energy) and F^^{x — )• 1) — > — oo while F^^{y q) ^ oo and F^^{y — )■ 1) — )■ —1. Note 
that at the two ends Fq{q) = —dq and -Fg(l) = oo. In spite of their difference, both 
F^^iz) and F^'^iz) have the same characteristic behavior as that of rriqiz), and this can be 
understood from the following relation (note its difference from ( fT4|) ). 



dF'^'^'^^z) ^ dmqjz) 
dz dz ' 

where we have used ( fT3|) . (fTS!) and the above Fq{z) (In the case of bubble, we also need 
to consider the contribution from dR/dy = dRq{y)/dy). For example, the maximum or 
minimum (if exists at all) for each of these functions occurs at the same ^max or 2min- The 
above enables us, in the present context, to make use of the extrema (if any) and behavior 
of mq{z) studied in [7j to compare the free energies at the respective global minimum. 

We are interested in the region for which the corresponding free energy will be at 
least locally stable, i.e., dmq{z)/dz\z=z < 0. Equation ( IT71) immediately tells us that 
the on-shell F^^{z) and F^^{z) along with mq{z) will decrease in this region. For rriqlz), 
b = bq{x) and R = Rq{y) then imply that in each of the cases b > R, b = R and b < R 
we will get x < y, x = y and x > y, respectively. Further, we can use the property of 
either F^^{x) or F^^{y) as a decreasing function to determine which of the phases - the 
brane or the bubble have smaller free energy for given b and R. Let us take a particular 
case for illustration, say b > R, implying x < y. So we have F^^{x) > F^^{y) = F^^{y), 
which says now that the bubble has smaller free energy. Alternatively, this can also be 
determined using F^*" as follows: F'°^{y) < F'^^(x) = -bFg{x) < -RFq{x) = F^'(x). In 
the following two sections, we will use F^^{z) as a decreasing function in the region of 
interest to determine the global stable phase, which appears slightly more straightforward 
than using F^^{z) instead. 

With the above preparation, we are now ready to discuss the phase structure for 
various cases as promised in the Introduction. 
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Figure 1: F^^{z) vs z and mQ{z) vs z for zero flux/charge case. 

3 The zero flux/charge case 

Let us discuss the zero charge/flux case first. As discussed in [7j, the function mo{z) 
at the two ends z = 0, 1 (look at eq. (fT3l) ) and has a maximum 



rrir. 



d 



in between at 2;max = 2/((i + 2) (see the second graph in Fig. 1). So, above mmax there is 
no black brane or bubble phase and the system will be in 'hot flat space' phase (with zero 
free energy). But below mmax, the fh = mo{z) always gives a local minimum of free energy 
at each Z2 lying in z-^i^^ < Z2 < 1 since there dmQ{z)/dz\z=z2 < 0, implying a locally stable 
black brane or bubble. The locally stable configuration becomes globally stable when 

(i.e, bg = Rg = nig) at Z2 lying in 1 > Z2 > Zg with 

. = (20) 

since its free energy is now negative. Here Zg is determined by setting the free energy to 
zero at Z2 = Zg > ^max- In other words, the free energy has its positive maximum value at 
^max, decreases from there to zero at Zg{> Zmax), then continues to decrease to negative 
and finally reaches — oo at ^2 = 1 (see the first graph for black brane in Fig. 1). So from 
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^2 = -^max = 2/((i + 2) to ^2 = 1, the free energy is a monotonically decreasing function 
and so, greater Z2 (i.e., X2 or 1/2) will give smaller free energy for the black brane or the 
bubble. Keeping this in mind we consider the following cases: 1) li R, b > bg = Rg, the 
black brane or bubble phase is at most locally stable with a positive free energy and the 
locally stable phase will make a Hawking-Page phase transition to the globally stable 'hot 
flat space' phase. If R > b = bg (or b > R = Rg), the bubble (or the black brane) is at 
most locally stable and the locally stable phase will make also a Hawking-Page transition 
to the globally stable phase which is now a coexisting phase of both the black brane (or 
the bubble) and the 'hot flat space'. 2) li R = b = bg = Rg, the three phases - the bubble, 
the black brane and the 'hot flat space' coexist and the transition between any two of 
them is a first order one since the first derivative of free energy has a discontinuity when 
the transition occurs. 3) If -R < Rg and R <b, the bubble is the globally stable phase. If 
b < bg and b < R, black brane is the globally stable phase instead. \ib = R<bg = Rg, 
the bubble and the black brane coexist and are globally stable. 



4 The non-zero flux/charge case 

We now move on to the case when the charge enclosed in the cavity is fixed but non-zero. 
As shown in [7], when the charge is non-zero, there exists a critical charge qc at which 
the first and second derivatives of mq{z) with respect to 2; at 2; = vanish and these will 
determine completely Qc, and rric. The underlying phase structure crucially depends on 
whether q > Qc, Q = Qc and q < Qc, so we discuss them in turns in the following for d > 2 
(The d < 2 cases are different and will be discussed afterwards.): i) For q > q^, there is 
no extrema for either mq{z) or free energy F^^'^^{z), with each decreasing monotonically 
in the region oi q < z, z < 1 (see the first graph in Fig. 2 for mq{z) and in Fig. 3 for 
F^^{z)). The slope of mq{z) is always negative and therefore the on-shell F^^'^^{z) is a 
local minimum. So, we have the following cases to consider. If 6 > ^, the respective 
rh = ruqiz) implies x < y and so, F^^{x) > F^^{y), i.e., bubble is the globally stable 
phase. In other words, an analog of Hawking-Page transition will take the locally stable 
black brane phase to the globally bubble phase. Similarly, i{b<R, the black brane is the 
globally stable phase. But iib = R, the bubble and the black brane phases coexist and the 
transition between the two is a first order one due to the change of entropy involved, ii) 
For q = qc, we have be = Rc and Xc = Vc, and the two phases are both globally stable and 
can coexist at the critical point and the transition between the two is also a first order one. 
iii) For q < q^, the situation is a bit involved. As noted in [7], in this case, the function 
mq{z) does not decrease monotonically in the region q < z < 1, but in between there is a 
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Figure 2: The typical behavior of mq[z) vs z for q > Qc and q < q, 



minimum m^n^ at z = Zram and a maximum mmax at 2 = z^^x (see the second graph in 
Fig. 2). Since mg{z) starts at infinity at 2; = g and goes to zero at 2; = 1 (see eq. (fT3|) ). 
so 2;min < -^max- Then as shown in [7], in the range q < z < z^ij^, and 2max < z < 1 the 
black branes or the bubbles are locally stable since each z gives a local minimum of free 
energy, whereas in z^i^i < z < ^max they are unstable since each z gives a maximum of free 
energy. So, for mmin < m < mmax, 'fn = mq{z) gives three black brane or bubble phases 
with three solutions, say, zi < ^2 < ^3, where zi (small black brane or bubble) and z^ 
(large black brane or bubble) correspond to the locally stable phases and 22 corresponds 
to unstable phase (see the second graph in Fig. 2). Among the locally stable phases the 
system would prefer to be in the phase of lowest free energy. Now as demonstrated in [7] , 
in between mmax and mmin and for given g, there exists a m denoted by mt{q,d) (here 
bf = Rt), a function of d and charge q only, where the large black brane (bubble) of size x^t 
(yst) has the same free energy {F^^'^^{zst) = F^^'^^{zit)) as the small black brane (bubble) 
of size Xit [yit) and so they coexist. Above rrit and below mmax the small black brane 
(bubble) is globally stable and below nit and above mmin, the large black brane (bubble) 
is the globally stable phase if the black brane or the bubble is assumed to be the only 
phase. Given what has been said about either black branes or bubbles, we now determine 
the globally stable phase with the lowest free energy between black branes and bubbles 
with the same boundary data. Given that F^^'^^{z) decreases monotonically in the range 
q < z < 2inin and ^max < z < 1 (see the second graph in Fig. 3), so when b, R > bt = Rt 
OT b, R < bt = Rt, which phase is globally stable can be discussed following q > qc case 
and will not be repeated here. When b = R = bt = Rt, all four phases i.e., the small black 
brane, the large black brane, the small bubble and the large bubble coexist. The transition 
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Figure 3: The typical behavior of reduced free energy F (z) vs 2 for g > and q < qc- 

from one phase to the other is a first order one (when it is between small and large black 
branes or between small and large bubbles) ending at a second order critical point when 
q = qc and a first order one (when it is between a black brane and a bubble) even at q = qc 
(due to a change of entropy there). The reason that this transition is always a first order 
one when q < qc is again due to the discontinuity of the first derivative of free energy 
when the transition occurs. When b > bf = Rt and R < bt = Rt, x < Xu < x^t = list < y 
and so, F^\x) > F^\xit) = F^\xst) = F'°^{y3t) > F'^^iV), therefore, the large bubble is 
the globally stable phase. Similarly, when b < bt = Rt and R > bt = Rt, the large black 
brane is the globally stable phase. 

The extremal brane phase with the same boundary data has always the largest free 
energy in any of the above three cases and therefore can never be a globally stable phase. 

We now move on to the d < 2 cases which have similarities with the zero charge case. 
Let us consider d = 1 case first. The phase analysis for this case can be discussed following 
the chargeless case if we replace the 'hot flat space' there by the present extremal brane 
phase and 1 > 2; > for mo (2;) there hj 1 > z > q for mq{z) here and will not be repeated 
here. The present mmax and Zmax can be determined now from mq{z) similarly but they 
cannot be given analytically as demonstrated in [7] . Instead, we will give their numerical 
values for a few selected values of q along with the values of rrig and Zg in the following 
table. The Zg falling in the range 1 > Zg > ^max can be determined from the following 
equation, 

which is obtained by equating the system free energy aX z = Zg with the corresponding 
extremal brane free energy at z = q, i.e., Fg{z = Zg) = Fq{z = q) = q- This is actually a 
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quartic equation but one can check, as expected, that 2; = g is a solution. So this equation 
can be reduced to a third-order equation as 

9z^ - {hq^ -q + 8)z^ + (3g + 4:)q'^z - V = 0, (22) 

which can be solved analytically and has only one real positive solution, giving Zg. The 
solution for Zg is complicated and not very illuminating and for this reason we will not 
give its explicit analytic expression here. As mentioned above, we will list instead its 
value for each given q along with the corresponding nig for q starting at g = 0.1 with an 
increment of 0.1 up to g = 0.9 in the table below. 



Q 




nig 




^raax 


0.1 


0.878471 


0.307757 


0.668371 


0.386344 


0.2 


0.870088 


0.319918 


0.673943 


0.390683 


0.3 


0.864945 


0.332758 


0.673943 


0.397903 


0.4 


0.864458 


0.346219 


0.684781 


0.407890 


0.5 


0.870015 


0.360200 


0.731065 


0.420323 


0.6 


0.882600 


0.374565 


0.769645 


0.434666 


0.7 


0.902505 


0.389161 


0.817842 


0.450290 


0.8 


0.929316 


0.403845 


0.873604 


0.466636 


0.9 


0.962160 


0.418491 


0.934861 


0.483295 



From the above table, one can see clearly that 1 > Zg > z^a.x and < mmax as expected. 

Let us now discuss d = 2 case. For this case, as shown also in [7], there exists not only 
(like d > 2 case) a critical charge withlj qc = 1/3, but also (like d = 1 case) a maximum 
of niq{z) at 2;max- For q > qc, have mmax = "^g(Q') = a/^/2 at z^a_^ = q (see the first 
graph in Fig. 4). When rfi > m^^^, the extremal brane is the global phase. But below 
mmax, this case is similar to the analysis oi d > 2 and will not be repeated here. When 
g < gc = 1/3, we have now 

1/2 



niy 



(l + 3g2 + V(l-g2)(l-9g2)) (^3 v^T^ - ^1^9?) 



8v/2 

1 + 3g2 + v/(l-g')(l-9g2) 



(23) 



^ However, now Zc = qc, occurring at the small end, and therefore, unlike d > 2, there is now no 
critical behavior since we don't have a stable small brane (or bubble) phase. 
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^ z 



Figure 4: The typical behavior of mq{z) vs z for g > Qc and q < qc when d = 2. 

with mjnax > iTT'qiq) = and Zmax > 5' (see the second graph in Fig. 4). This subcase 

is however similar to the discussion of d = 1. Again when rh > m^^i^^, the extremal brane 
is the stable phase. When b < b^ax and R > -Rmax (or b > femax and R < -Rmax) there exists 
a locally stable black brane (or bubble). As in (i = 1 case, we now need to determine 
Xg (or yg) below which locally stable black brane (or bubble) becomes globally stable. 
Exactly as in d = 1 case, Zg is now determined by the following equation, 

— (Tfi)""*(a)"*.-'"(-f)" 

This equation, unlike in the previous case, can be solved exactly and has apart from the 
trivial solution at ^ = g, two other real solutions. Only one of them lies in the range 
2^max < < ^ and is given as 

_ 3 - 2g + + (1 + g) ^3(3 - g)(l - 3g) 

with 1 > Zg > Zmax With this Zg, we have 

5 - 3g - v/3(3-g)(l -3g)j [3 - 2g + ?>q^ + (1 + g) ^3(3 - g)(l - 3g)j ^'^ 



rUg 



16\/2 

(26) 

Now the discussion for globally stable phase and the analog of Hawking-Page transition 
are identical to the chargeless case if we replace the 'hot flat space' there by the present 
extremal brane phase and therefore will not be repeated. 
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5 Discussion and conclusion 



In this paper, we find that a necessary completion of phase structure of charged black 
p-branes in canonical ensemble requires two additional phases, namely, the bubble and 
the extremal brane. This finding solves the puzzle about the missing phase(s) for d < 2 
and gives a new phase for d > 2. We also find an analog of Hawking-Page transition 
among many new phase transitions revealed, giving now a much enriched and complete 
underlying phase structure (assuming no other new phases present) which can have a 
coexistence of up to four individual phases. The phase transitions are always the first- 
order ones. We also find that the extremal brane cannot be a stable phase for d > 2 but 
is vital to the completion of phase diagram for < 2. These results are obtained solely 
on the basis of the free energy criterion which is used to justify the underlying stability 
of the various phases. 

Before closing this section, we need to clarify which spin structure the extremal brane 
should take so that no inconsistency arises. For this, let us be somewhat specific in the 
present context. We focus on the relevant {t,x^,p) directions. For the non-extremal 
branes, the topology is x with the denoting the x^-circle while for the bubble it 
is 5*^ X i?^ with now denoting the t-circle. Though the two have the same topology in 
Euclidean signature, they are quite different in Lorentzian signature since their spacetime 
structures are different. Also the former has a non- vanishing entropy while the latter 
has zero entropy. In spite of their differences, the two share one common feature that 
there is only one spin-structure allowed corresponding to the fermions being anti-periodic 
along their respective S^. However, for the extremal brane, the topology is x x R 
and in general there are two spin structures, one for fermions being periodic and the 
other for fermions being anti-periodic, along either circle. The extremal p-branes, when 
considered in isolation, are the 1/2-BPS objects in string/M theory, which preserve one 
half of the spacetime super symmetries. So the natural choice for the spin structure is with 
periodic fermions. If this were the case in the present context, one would have difficulty in 
understanding^ the phase transitions between the extremal brane and the non-extremal 
brane (or bubble) for d < 2, vital for the completion of phase diagram, since the two have 
different spin structure as discussed above. However, the present setup of placing the 
extremal brane in a cavity with a fixed finite temperature in canonical ensemble breaks 
all the underlying supersymmetries and therefore requires the spin structure with anti- 

^No such issues would arise for d > 2 since there is no phase transition between the extremal brane 
and the non-extremal brane (or the bubble). 
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periodic fermions [20] for the extremal branqj instead. So spin structure wise, there is no 
paradox arising and the free energy criterion is as usual a suitable means to obtain the 
phase structure and phase transitions as given in this paper. 
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